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Abstract 

The X(3872) seems to be a loosely-bound hadronic molecule whose constituents are two charm 
mesons. A novel feature of this molecule is that the mass difference of the constituents is close 
to the mass of a lighter meson that can be exchanged between them, namely the pion. We ana- 
lyze this feature in a simple model with spin-0 mesons only. Various observables are calculated to 
next-to-leading order in the interaction strength of the exchanged meson. Renormalization requires 
summing a geometric series of next-to- leading order corrections. The dependence of observables on 
the ultraviolet cutoff can be removed by renormalizations of the mass of the heaviest meson, the 
coupling constant for the contact interaction between the heavy mesons, and short-distance coeffi- 
cients in the operator product expansion. The next-to-leading order correction has an unphysical 
infrared divergence at the threshold of the two heavier mesons that can be eliminated by a further 
resummation that takes into account the nonzero width of the heaviest meson. 

PACS numbers: 12.38.-t, 12.39. St, 13.20. Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a new hadronic resonance discovered in 2003 by the Belle Collaboration 
[1] and subsequently confirmed by the CDF, Babar, and DO Collaborations [2-4]. All the 
properties of the X(3872) that have been measured thus far are compatible with its iden- 
tification as a weakly-bound molecule whose constituents are a superposition of the charm 
meson pairs D*°D° and D°D*° [5-25]. 

Because it is so weakly bound, the X(3872) has many features in common with the 
deuteron, which is a weakly-bound baryonic molecule consisting of a proton and a neutron. 
Their binding energies are both small compared to the natural energy scale associated with 
pion exchange: m^/(2Mi2), where M\ 2 is the reduced mass of the two constituents. The 
binding energy 2.2 MeV of the deuteron is small compared to the natural scale of about 
20 MeV. After taking into account a recent precision measurement of the D° mass by the 
CLEO Collaboration [26], the difference between the measured mass of the X(3872) and the 
D*°D° threshold is 

M x - (M 1 + M 2 ) = -0.6 ± 0.6 MeV. (1) 

This is small compared to the natural energy scale of about 10 MeV. The Belle Collaboration 
has set an upper bound on the width of the X [1]: 

Y x < 2.3 MeV (90% C.L.). (2) 

This is also small compared to the natural energy scale. The small width is most easily 
understood if the X is below the D*°D° threshold. 

Because the binding energies of the deuteron and the X(3872) are small compared to 
the natural energy scales, they have universal properties that are determined by the large 
scattering length a of the constituents [10]. The universality of few-body systems with a 
large scattering length has many applications in atomic, nuclear, and particle physics [27]. 
The universal features of the X(3872) were first exploited by Voloshin to describe its decays 
into D°D 7i° and D°D 0/ -f, which can proceed through decay of the constituent D*° or D*° 
[8]. Universality has also been applied to the production process B — > KX [14, 17], to the 
line shape of the X [18], and to decays of X into J/ip and pions [20]. These applications 
rely on factorization formulas that separate the length scale a from all the shorter distance 
scales of QCD [18]. The factorization formulas can be derived using the operator product 
expansion for a low-energy effective field theory [24]. 

There is one feature of the X(3872) that is very different from the deuteron. The mass 
difference between the constituents D*° and D°, M 2 -M 1 = 142.12±0.07 MeV, is very close 
to the 7T° mass: m = 134.98 MeV. Thus the D*° mass is very close to the D n threshold: 
M 2 — Mi —m = 7.14 ± 0.07 MeV. This is important because the pion is the lightest meson 
that can be exchanged between the D*° and D°. The small energy gap between the D*° 
mass and the D°7r° threshold makes it easy for the D*° to make a transition to D 7t°. An 
important consequence is that the X could have a substantial 3-body component D°D°7r°. 
Suzuki has argued that the near equality of M 2 — Mi and m implies that the pion-exchange 
interaction is too weak to bind the X [21], but this conclusion has been criticized [23]. 

An analysis of the effect of the near equality M 2 — M 1 ps m in the case of the X(3872) is 
complicated by other effects that may be comparable in importance. There are other nearby 
two-body and three-body thresholds. The thresholds for the charged mesons D* + D~ and 
D+D*- are higher than the D*°D° threshold by 8.08±0.12 MeV. The D°D'n + and D + D°u~ 
thresholds are higher than the D*°D° threshold by only 2.23 ± 0.12 MeV. There are several 



2 



other features that complicate the analysis of the effects of pions on the X. The constituent 
D*'s are spin-1 particles. The charge conjugation C = + requires the molecule to be a 
superposition of charm mesons: D*°D° + D°D*°. The interaction of the charm mesons with 
pions depends on the pion momentum. A quantitative analysis of the effects of pions on the 
X might require all these effects to be taken into account. 

Since the new feature of the X(3872) is the near equality between the mass difference 
M 2 — Mi of the constituents and the mass m of a meson that can be exchanged between 
them, it is worthwhile to analyze this feature in a simple model that avoids all the other 
complications of the X. The simplest such model is one in which the primary constituents 
D\ and D 2 are spin-0 mesons with a large positive scattering length and a momentum- 
independent coupling to a lighter spin-0 meson <fi. We will calculate the effects of the 
exchanged meson in this model to next-to-leading order in the D 2 — Di<p coupling constant 

9- 

In Sec. II, we write down the lagrangian for the model and discuss its parameters. In 
Sec. Ill, we summarize the results of the model at th order in g that were derived previously 
in Refs. [18, 24]. The results include the short-distance production rate for X and the 
line shape of X in a short- distance decay channel. The ultraviolet divergences can be 
removed by renormalization of the coupling constant for the D\D 2 contact interaction and by 
renormalization of short-distance coefficients in the operator product expansion. In Sec. IV, 
we calculate the observables in Sec. Ill to next-to-leading order in g. We also calculate the 
short-distance production rate for D 1 D 1 (f) and the decay rate of X into DiDi<p. There are 
new ultraviolet divergences that can be removed by renormalization of the D 2 mass. We 
show that after summing a geometric series of next-to-leading order corrections, all remaining 
ultraviolet divergences can again be removed by renormalization of the coupling constant 
for the D\D 2 contact interaction and by renormalization of short-distance coefficients in the 
operator product expansion. A further resummation is required to eliminate an infrared 
divergence at the D\D 2 threshold in the next-to- leading order corrections. In Sec. V, we 
summarize our results. The two-loop Feynman diagrams that arise at next-to-leading order 
in g are calculated in an Appendix. 

II. A SCALAR MESON MODEL 
A. Lagrangian 

We consider the simplest possible model with a bound state X whose mass is close to a 
3-body threshold. The primary constituents D\ and D 2 of the bound state are scalar mesons 
whose masses M x and M 2 satisfy M x < M 2 . There is also a scalar meson whose mass m is 
close to the mass difference M 2 — M x . We will refer to this model as the scalar meson model. 
It is defined by a nonrelativistic quantum field theory with three complex scalar fields D\, 
D 2 , and 0. The free terms in the lagrangian 1 are 

i=l,2 x 7 x ' 



The superscript f on a field indicates its complex conjugate. 
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FIG. 1: Vertices for (a) the D\D 2 contact interaction, (b) the D\D 2 (p interaction, and (c) the D 2 
mass counterterm. 

The interaction terms for the scalar meson model are 

C mt = -XoDlDlD^ - g(DlD 1( p + D\^D 2 ) - 8MD\d 2 . (4) 

The interaction vertices are illustrated in Fig. 1. The coupling constants Ao and g have mass 
dimensions —2 and —1/2, respectively. We assume that there is a fine-tuning that makes the 
-D1-D2 scattering length large, which implies that the D\D 2 contact interaction with coupling 
constant Ao must be treated nonperturbatively. The subscript on Ao emphasizes that it is 
a bare coupling constant that requires renormalization. We assume that the D 2 — Di<j) 
interaction can be treated perturbatively and we calculate its effects through order g 2 . At 
this order, no renormalization of g is required. However, we will see that nonperturbative 
renormalization of Ao requires summing a geometric series to all orders in g 2 . There are 
also ultraviolet divergent corrections to the mass of the D 2 which are cancelled by the mass 
counterterm 5M in Eq. (4). 



B. Masses 

It is convenient to introduce concise notations for various combinations of the masses. 
The sum and difference of the D\ and D 2 masses are 

M 1+2 = M 1 + M 2l (5a) 
M 2 _i = M 2 -M 1 . (5b) 

The reduced masses of D\D 2 and of D\(j) are 

12 M l + M 2 y 1 

M l + m y ' 

If m is much smaller than Mi, the natural energy and 3- momentum scales associated with 
exchange of the meson between the mesons D\ and D 2 are m 2 /M\ 2 and m, respectively. We 
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assume that the energy gap between the D 2 mass and the D\<\> threshold is small compared 
to the natural energy scale: 

M 2 _i -m < m 2 /M 12 . (7) 

The total number of heavy mesons D 1 and D 2 is conserved by the interactions in Eq. (4). 
We are interested in the sector of the theory with two heavy mesons and, more specifically, 
in the threshold region where the invariant mass M of all the particles is very close to M 1+2 : 

\M-M 1+2 \ <m 2 /M 12 . (8) 

This restricts the possible scattering states to D X D 2 and DiDi<p. 

For numerical illustrations, we will use masses that correspond to the D*° D° / D° D°7r° 
system. The meson masses are M 2 = 2006.7 MeV, Mi = 1864.5 MeV, and m = 135.0 MeV. 
Thus M12 « 966 MeV and M 2 _i - m = 7.2 MeV. The natural energy scale is m 2 /M 12 « 19 
MeV. The condition in Eq. (7) is only marginally satisfied: Mi 2 (M 2 _i — m)/m 2 0.4. We 
take the binding energy of X to be Ex = 0.4 MeV, which is near the central value of the 
measurement of Mx in Eq. (1). If we express this binding energy as Ex = k x /(2Mi 2 ), the 
binding momentum is Kx = 27.8 MeV. 



C. Coupling constants 

We impose an ultraviolet cutoff on the momenta p of particles in the center-of-momentum 
frame that restricts them to the region \p\ <C m. This guarantees that all the particles 
remain nonrelativistic. It was shown in Ref. [24] that at leading order g°, all ultraviolet 
divergences can be absorbed into A and into short-distance coefficients in the operator 
product expansion. At order g 2 , there are new ultraviolet divergences that can be absorbed 
into the D 2 mass. We will find that the remaining ultraviolet divergences can again be 
absorbed into Ao and into short-distance coefficients in the operator product expansion. 

The coupling constant Ao must be tuned to near the critical value at which the D\D 2 
scattering length diverges in order for the bound state X to have a small binding energy 
Ex <C m 2 /M\ 2 . In order to compare results at order g° and order g 2 , we must have a well- 
defined renormalization prescription for Ao or, equivalently, a clear definition of the binding 
energy Ex- We will find that all the amplitudes for processes near the D\D 2 threshold 
have a resonant term with a common factor A(E) that depends on the energy E in the rest 
frame. For example, the line shape in a short-distance decay channel of X produced by a 
short-distance process has a long-distance factor |„4.(M)| 2 , where M is the invariant mass of 
the decay products. The square of the resonant factor can be expressed in the form 

|,A(jB)|2 = [ReA' l (E)] 2 + [Im^- 1 ^)] 2 ' (9) 

As a function of the real energy E, this resonance factor has a peak when E is close to 
Mi +2 — Ex- The peak arises because Rev4. _1 (-E) vanishes near that value of E and lmA~ 1 (E) 
is much smaller than Re^4 _1 (-E) except near that value of E. As our renormalization 
prescription for A , we demand that Re^4~ 1 (£') vanishes at an energy determined by a real 
variable k X '- 

Re^- 1 ^) = at E = M l+2 - k 2 x /{2M 12 ). (10) 
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The dependence of an observable on the bare parameter Ao can be eliminated in favor of 
the renormalization parameter Kx- We assume that Kx is small compared to the natural 
momentum scale: Kx <C m. 

The amplitude A(E) has a pole at a complex energy -E po i e near the real energy at which 
|„4(i?)| 2 has its maximum. We choose to define the mass M x and the width T x of X by 
expressing that complex energy in the form 

£ P oic = M x - iTx/2. (11) 

We also choose to define the binding energy E x by 

E x = M 1+2 - M x . (12) 

We assume that the binding energy and the width of X are both small compared to the 
natural energy scale: E x , ^x <C m 2 /M 12 . 

In Ref. [24], the authors proposed a renormalization scheme in which Ao was tuned to 
give a prescribed value for the complex parameter E po \ e . The binding energy Ex and the 
width Tx were taken as the input parameters that determine E po \ e . This renormalization 
scheme requires Ao to be a complex coupling constant. The imaginary part of Ao takes into 
account inelastic scattering channels for DiD 2 that are outside the energy range described 
by the effective field theory. In this paper, we assume that any such inelastic channels do 
not give a significant contribution to the width of X. The coupling constant Ao is therefore 
real valued. We could use the binding energy E x as the input parameter and eliminate A 
in favor of E x , and then T x would be calculable. We find it more convenient to eliminate 
Ao in favor of the variable Kx introduced in Eq. (10). In this renormalization scheme, Ex 
and Tx are both calculable in terms of kx, 9, and the masses Mi, M 2 , and m. 

The most convenient input for determining the coupling constant g is the width r 2 of the 
heavy meson D 2 . Since M 2 > Mi + m, the DiD 2 (f) vertex in Fig. 1(b) allows D 2 to decay 
into Di4>. Calculating the width using nonrelativistic phase space, we obtain 

r 2 = ^i[ 2mi (M 2 _ 1 -m)] 1/2 , (13) 

7T 

where Mi 2 and mi are the reduced masses defined in Eqs. (6). The value of the coupling 
constant g can be determined from T 2 and the masses. We assume that the width of D 2 
is small compared to the natural energy scale: T 2 <C m 2 /M 12 . We also assume T 2 is small 
enough to allow perturbation theory in the coupling constant g. 

For numerical illustrations, we will use values of the parameters that correspond to the 
D*°D /D D 7i° system. We take the width of D 2 to be Y 2 = 0.07 MeV, which is ap- 
proximately the total width for D*° [18]. This is much smaller than the natural energy 
scale m 2 /M 12 and also rather small compared to M 2 _i — m: T 2 /(M 2 _! — to) ~ 0.01. Using 
Eq. (13), we determine the numerical value of the DiD 2 (p coupling constant to be g = 0.0064 
MeV" 1 / 2 . When considering the dependence of observables on g, it will be more convenient 
to regard them as functions of T 2 . 

We now list all the energies that are assumed to be small compared to the natural scale 
m 2 /M 12 : 

• the energy M — M 1+2 of a D X D 2 or DiDi<j) system relative to the DiD 2 threshold, 

• the energy difference M 2 _i — to between the D 2 mass and the Di4> threshold, 
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• the decay width T 2 of D 2 , which is given in Eq. (13), 

• the binding energy E x of the molecule, which is approximately equal to k|/(2M 12 ), 
where kx is a renormalization parameter. 

We will carry out our calculations under the assumption that these energies are all compa- 
rable and that the masses M 2 , Mi, and m are also comparable. The only approximations 
we will make are those that can be justified by a nonrelativistic approximation, such as 
M 2 _x - m < M 12 . 



D. Self-energy of 



At order g 2 , renormalization of the D 2 mass is necessary. In the on-shell renormalization 
scheme, the parameter M 2 in the free lagrangian in Eq. (3) is chosen to be the physical D 2 
mass. The physical mass M 2 and the width T 2 of D 2 can be defined by specifying that the 
exact D 2 propagator at momentum has a pole in the energy at po = M 2 — iT 2 /2. The pole 
in the energy po at momentum is the solution to the equation 



Po -M 2 - S(p , 0) = 0, 



(14) 



where S(p ,p) is the D 2 self-energy. The exact D 2 self-energy is the sum of the one-loop 
diagram in Fig. 2(a) and the mass counterterm in Fig. 2(b): 



£(Po,p) = £ 2 a(po,p) + SM. 



(15) 



There are no other diagrams for the self-energy at higher order in g. The one-loop diagram 
is calculated in Section A 2 of the Appendix. The analytic expression for the self-energy 
diagram using dimensional regularization is given in Eq. (A9). The result in a general 
regularization scheme is 



£ 2 a(po,p) = S 2a (Mi +m, 0) -% 



.g 2 m 1 



2n 



2mi 



Po — Mi — m — 



p 



2(Mi + m) 



+ ie 



1/2 



(16) 



The term S 2a (Mx + m, 0) is real valued and includes a linear ultraviolet divergence. The 
divergence is cancelled by a divergent term in the mass counterterm 5M in Eq. (15). In 
the on-shell renormalization scheme, the finite terms in SM are chosen so that the solution 
to Eq. (14) satisfies Repo — M 2 . The wavefunction normalization factor for a D 2 with 
momentum p is 

Z^p)' 1 = 1 - ^-S(M 2 - *r 2 /2 + p 2 /(2M 2 ),p). (17) 
opo 

Since Eq. (14) cannot be solved analytically, the on-shell renormalization scheme is not 
the most convenient choice. A simpler renormalization scheme for the D 2 mass is to choose 
the mass counterterm so the self-energy S 2a (p 0) 0) vanishes at the DiD 2 threshold: SM = 
— S 2a (M! + m, 0). The self-energy is then 



£(po,p) = 



.g 2 mi 



2mi 



Po - Mi - m 



p* 



2(Mi + m) 



+ ie 



1/2 



(18) 
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FIG. 2: The diagrams for the Z?2 self-energy. 

With this renormalization prescription, the parameter M 2 differs at order g from the phys- 
ical D 2 mass, which is given by the real part of the solution to Eq. (14). The width r 2 , 
which is given by the imaginary part of the solution to Eq. (14), reduces at leading order in 
g to T 2 = — 2 Im H 2a (M 2 , 0). This agrees with the explicit expression for the decay rate of 
D 2 at order g 2 obtained in Eq. (13). To order g 2 , the wavefunction normalization factor in 
Eq. (17) is 

Z 2 ( P y l = l + t^—^(2m 1 (M 2 _ 1 -m)[l-p 2 /(2(M 1 + m)M 2 )]) . (19) 

Since the momentum p is restricted to the nonrelativistic region, the term p 2 / (2(Mi +m)M 2 ) 
can always be neglected compared to 1. The resulting wavefunction normalization factor is 
independent of p: 

2 2 

Z- 1 = 1 + i 9 -^ [2m 1 {M 2 . 1 - m)]- 1/2 . (20) 
The correction term is pure imaginary. 



E. The Complex Mass Scheme 

The D 2 — D\(j) interaction allows a virtual D 2 whose energy is greater than M\ + m to 
decay into Di<fi. In the amplitude for such a process at leading order in g, the propagator 
of the virtual D 2 diverges. Near such a point in momentum space, the corrections to the 
D 2 propagator from higher orders in g are not suppressed. This problem can be solved by 
summing the D 2 propagator corrections to all orders. The effect is to add — H(p ,p) to the 
denominator po —p 2 / (2M 2 ) +ie of the D 2 propagator. The imaginary part of the self-energy 
S(po,p) eliminates any divergence in the D 2 propagator for real values of the energy. One 
drawback of this solution is that it leads to unnecessarily complicated expressions for the 
integrals in loop diagrams with D 2 propagators. 

Summing the D 2 self-energy corrections to all orders is essential only near the pole in 
the propagator. In other regions of the energy p Q , S(po,p) is a perturbative correction to 
the denominator of the propagator that is suppressed by a power of g 2 . In those regions, 
summing the D 2 propagator corrections to all orders is optional. At the pole in p , the 
self-energy ~E(p ,p) reduces to — iT 2 /2. An alternative resummation that eliminates the 
divergence in the propagator for real values of the energy without changing the order in 
g of the truncation error is to sum only the — iT 2 /2 term in the self-energy to all orders 
and treat the remainder of the self-energy as a perturbation. An advantage of this partial 
resummation is that it leads to simpler expressions for the integrals in loop diagrams with 
D 2 propagators. 



A systematic method for implementing this partial resummation is the complex mass 
scheme [28]. This method has been used in Standard Model calculations at next-to- leading 
order [29-31]. The application of this method to the scalar meson model involves adding 
cancelling terms to the free and interaction terms in the Lagrangian: 

A£ frcc = (iT 2 /2)DlD 2 + (Z 2 l - l)D\ (i^ -M 2 + iF 2 /2 + ^-V 2 ) D 2 , (21a) 

A£ int = -{iY 2 /2)D\D 2 - (Z, 1 - l)D\ (i^ - M 2 + iY 2 /2 + ^V 2 ) D 2 , (21b) 

where T 2 is the D 2 width and Z 2 is the D 2 wavefunction normalization factor. To order g 2 , T 2 
and Z 2 are given by Eqs. (13) and (20). Since iT 2 /2 is pure imaginary and Z 2 is complex, the 
free Lagrangian consisting of the sum of Eqs. (3) and (21a) corresponds to a nonhermitian 
Hamiltonian. The effects of the nonhermiticity are cancelled exactly if the interaction terms 
in Eq. (21b) are calculated to all orders. At any finite order in perturbation theory, the 
complex mass scheme gives amplitudes that are uniformly accurate to the appropriate order 
in g, even if there is a virtual D 2 that decays into Di<p. The effect of adding Eq. (21a) to 
the free Lagrangian in Eq. (3) is to change the D 2 propagator: 

(22) 



p -M 2 - p 2 /(2M 2 ) + ie po- (M 2 - *r 2 /2) - p 2 /(2M 2 ) ' 

The effect of adding Eq. (21b) to the interaction Lagrangian in Eq. (4) is to change the 
counterterm vertex in Fig. 1: 



-i5M 



- -i(5M + *r 2 /2) - i{Z 2 l - 1) [p -M 2 + iT 2 /2 - p 2 /{2M 2 )] . (23) 

It is convenient to also rescale the fields D 2 (x) and D 2 (x) by the same complex factor Z^ 2 . 
This eliminates the factor of Z 2 from the numerator of the D 2 propagator in Eq. (22). It 

1 /2 

also multiplies all the interaction terms in Eqs. (4) and (21b) by Z 2 or Z 2 ' . Since we will 
use nonperturbative renormalization for the coupling constant Z 2 \q for the D\D 2 contact 
interaction, the factor of Z 2 only affects the value of the unphysical bare coupling constant 
Ao- Since Z 2 = l + 0(g 2 ), the other factors of Z 2 or Z^ 2 in the interaction vertices contribute 
only at order g 4 and higher. We will be calculating only to order g 2 , so we will ignore these 
factors of Z 2 and Z X 2 2 . 

The perturbation expansion in the complex mass scheme corresponds to ordinary per- 
turbation theory in g followed by a resummation to all orders of a subset of terms that are 
higher order in g. We will refer to calculations to order g° in this perturbation expansion 
as leading order (LO) in the complex mass scheme. We will refer to calculations through 
order g 2 in this perturbation expansion as next-to-leading order (NLO) in the complex mass 
scheme. 



F. Propagator for X 

We can use the scalar meson model to describe processes with asymptotic states that 
include not only the particles Di and 0, but also the bound state X. The local composite 
operator D\D 2 (x) has a nonzero amplitude to create X from the vacuum. Thus DiD 2 (x) 



9 



can be used as an interpolating field for X. It is more convenient to use the operator 
\qDiD2(x), because this is a renormalized operator whose matrix elements do not depend 
on the ultraviolet cutoff [24]. The corresponding propagator for X is 

A X (E,P) = J d^e^mXoDiD^XoDlDliOm. (24) 

— * — * 

The momentum 4- vector in the Fourier transform is P M = (E,P), where E and P are the 
energy and momentum of X. At P = 0, this propagator has a pole at the complex energy 
E po \ e . Near the pole in the energy, the behavior of the propagator at P = is 

A x (E,0)^ f X (25) 

The energy E po \ c determines the binding energy Ex and the width Fx of X through Eq. (11). 
The residue of the pole defines the wavefunction normalization factor Z x , which is complex 
valued. Because the propagator A X (E,P) has mass dimension —2, Zx has mass dimension 
-1. 

Strictly speaking, since X has a nonzero width T x , there are no T-matrix elements for 
processes involving X in the initial or final state, because the decay of the X prevents 
it from being a truly asymptotic state. However if Tx is sufficiently small, the X can 
propagate over a long time interval before decaying, and it can therefore be treated as a 
quasi-asymptotic state. We can use the Lehmann-Symanzik-Zimmermann (LSZ) formalism 
to define T-matrix elements for the X that are closely related to T-matrix elements for 
decay products of the X whose total energy is tuned to the peak of the X resonance. To 
define the T-matrix element T for a process with X in the initial or final state, we start 
with a Green's functions for the operator \qD\d\(x) or \ D 1 D 2 (x). The connected Green's 

— * 

function at 4-momentum (E, P ) is amputated by multiplying by the inverse propagator 
A X (E, P)" 1 , it is normalized by multiplying by Z][ 2 , and then it is evaluated at the real 
energy Mx + P 2 / (2M X ) ■ This gives the T-matrix element multiplied by % for a state X with 
the standard nonrelativistic normalization. The T-matrix element T for a state X with the 
standard relativistic normalization is obtained by multiplying by ^/2Mx- 

To justify this prescription, we consider a process A — > B + C, where A and B both 
represent one or more particles and C denotes a set of particles that can be decay products 
of X. The T-matrix for this process will have a resonant enhancement when the invariant 
mass M of the particles in C is near the mass Mx- The resonant contribution to the 
T-matrix element for A — > B + C in the rest frame of C can be approximated by 

iT[A ^B + C],iT [ A^B + X] — - ^ - ^ iT[X C]. (26) 

Our prescription for T-matrix elements such as T[A — > B + X] and T[X — > C] involves 
evaluating the amputated connected Green's function at the real energy E = Mx- This 
ensures that the expression in Eq. (26) is as good an approximation as possible to the 
T-matrix element T[A — > B + C] near the peak of the resonance. 

As a simple illustration of the prescription for T-matrix elements for processes involving 
X, we consider the forward scattering process X — > X. The relevant Green's function in the 
rest frame of X is the negative of the propagator A X (E, 0). This must be multiplied by two 
factors of A X (E, 0)~ 1 Z X /2 , one for the X in the inital state and one for the X in the final 
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state. The T-matrix element iT[X — > X] is then obtained by evaluating this at E = Mx 
and multiplying it by 2Mx- 



G. Short-distance production and decay processes 

The scalar meson model defined by the lagrangian in Eqs. (3) and (4) can be a low-energy 
approximation to a more fundamental quantum field theory. The fundamental theory may 
include high energy processes that can create DiD 2 or DiDi<f) with invariant mass near 
M 1+2 . If long-distance effects involving momenta much smaller than m can be separated from 
short- distance effects involving momenta of order m and larger, we can use the scalar meson 
model to calculate the long-distance effects. The basic tool required to separate long-distance 
effects from short-distance effects is the operator product expansion [24]. Applications of the 
operator product expansion to the scalar meson model at th order in g were described 
in Ref. [24]. Below, we give the leading terms in the operator product expansions of the 
T-matrix elements for some of the processes considered in Ref. [24]. 

The general production process for X has the form A — > B + X, where A and B each 
represent one or more particles. There can be analogous production processes for D\D\§. In 
the case of DiDi<p, we assume that the invariant mass M satisfies the condition in Eq. (8). 
This implies that the relative momenta of the mesons, which we denote generically by p, 
are small compared to the natural momentum scale m. We call the production process 
a short- distance process if all the particles in A and B have momenta in the rest frame 
of X or D 1 D i (f) that are of order m or larger. The T-matrix element for a short- distance 
production process can be expanded in powers of the small energy differences M — M 1+2 and 
M 2 _i — m divided by m 2 /M 12 and higher energy scales and in powers of the small relative 
momenta p divided by m and higher momentum scales. The operator product expansion 
can be used to organize the expansions of the T-matrix elements into sums of products of 
short-distance coefficients and matrix elements of local operators between the vacuum state 
|0) and the appropriate final state. The leading terms in the expansions are those with the 
lowest dimension operator, which is D\d\(x). The corresponding renormalized operator is 
\qD\d\(x). In Feynman diagrams, the local operator \ D\d\(x) is represented by a dot 
from which a Di line and a D 2 line emerge, as illustrated in Fig. 3. If we keep only the 
leading terms, the factorization formulas for the T-matrix elements reduce to 2 

T[A^B + X] = y/2Mx~C*' 12 (X\X d\dI(0)\$) : (28a) 
T[A -> B + £>i£>i</>] = ^8M 2 m C B / 2 {D x D x $\\D\d\{$)\§) . (28b) 

The Wilson coefficient C^' 12 is a function of high energy scales, such as m and M i2 , and of 
the total momentum of the X or DiDi<p. The only dependence on whether the final state 
contains X or DiD^ is in the matrix elements of the local operator X o D[Dl(0). Short- 
distance and long-distance effects are separated in Eqs. (28). 



2 The Wilson coefficients C A ' 12 and differ from those in Ref. [24] by a factor of Ao- This choice simplifies 
many equations. 
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FIG. 3: Vertex for the XqD\D 2 and \qD\D\ operators. 

The decay of X into D\D\(j) can be described within the scalar meson model. The 
fundamental theory may also allow the decay of X into other final states C. We call such 
a process a short-distance decay if all the particles in C have momenta in the rest frame 
that are of order m or larger. The T-matrix element for such a process can be expanded 
in powers of the small energy differences Mx — M\ +2 and M 2 -x — m divided by m 2 jM\ 2 
and higher energy scales. The operator product expansion can be used to organize the 
expansion of the T-matrix element into sums of products of short-distance coefficients and 
matrix elements of local operators between the state \X) and the vacuum state. The leading 
term in the expansion is the one with the lowest dimension operator, which is DiD 2 (x). 
The corresponding renormalized operator is XoDiD 2 (x). If we keep only this term, the 
factorization formula for the T-matrix element reduces to 2 

T[X C] = y/2M^Cg (®\\ D 1 D 2 (0)\X). (29) 

The Wilson coefficient Cf 2 is a function of high energy scales, such as m and M\ 2 . Short- 
distance and long-distance effects are separated in Eqs. (29). 

If the fundamental theory includes processes that allow the production of X via A — > 
B + X and the decay of X via X — > C, it also allows the process A — ► B + C, where C 
represents the same particles but with a variable invariant mass M instead of Mx- This 
process has a resonant enhancement when M is near the D\D 2 threshold as specified by 
Eq. (8). We call it a short- distance process if each of the particles in A and B has momentum 
large compared to m in the rest frame of C and if the relative momentum between each 
particle in C and each particle in A or £ is large compared to m. The T-matrix element 
for such a process can be described within the scalar meson model by a double operator 
product expansion. The leading terms in this expansion are 

T[A^B + C] = Cf c + Cf' 12 CgJ d^e^mXoD.D^XoDlDliom, (30) 

where the 4- vector is P M = (M, 0). The Wilson coefficients C^' 12 and Cg are the same 
ones that appear in the operator product expansions in Eqs. (28) and (29). The first term 
C^' c on the right side of Eq. (30) takes into account the direct production of C at short 
distances. This term can be expanded in powers of the small energy differences M — M 1+2 
and M 2 _i — m divided by energy scales of order m 2 /M 12 or higher. The leading term in this 
expansion is simply a constant. According to Eq. (24), the Fourier transform in Eq. (30) is 
just the X propagator evaluated at E — M . Short-distance and long-distance effects are not 
yet separated in Eq. (30), because the X propagator requires an additive renormalization 
[24]. 
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FIG. 4: Diagram of order g° for the amplitude for the propagation of D\D 2 between contact 
interactions. The open dots indicate that the vertex factors — zAo are omitted. 

III. LEADING ORDER IN g 



At th order in g, the only interaction in the scalar meson model is the contact interaction 
between D\ and D 2 with coupling constant Ao- The DiDi<p states remain noninteracting at 
this order in g. In this section, we summarize some of the results at th order in g that were 
obtained in Ref. [24] and we generalize those results to the complex mass scheme. 



A. Amplitude for D^D? — > D^D^ 

At th order in g, all the observables for processes near the D\D 2 threshold are related 
in a simple way to the connected Green's function for D\D 2 — > D\D 2 . The basic building 
block for this Green's function is the amplitude iLq(E) for the propagation of the D\D 2 pair 
between successive contact interactions, which is given by the Feynman diagram in Fig. 4. 
The amplitude L (E) is calculated in Section A3 of the Appendix. The analytic expression 
using dimensional regularization is given in Eq. (A13). The result in a general regularization 
scheme is 

L (E) = L (M 1+2 ) + ^k(E), (31) 
where k{E) is the energy variable defined by 



k{E) = V-2M 12 (£ - M 1+2 ) - is. (32) 

This variable vanishes at the D\D 2 threshold E = M\ +2 . It is real and positive if E < M\ +2 
and it is pure imaginary with a negative imaginary part if E > M\ +2 . The term Lq(Mi +2 ) 
in Eq. (31) is real valued and includes a linear ultraviolet divergence. 

The amputated connected Green's function %Aq{E) for D 1 D 2 — > D X D 2 can be calculated 
nonperturbatively by summing the geometric series represented by Fig. 5 to all orders in Ao: 

^ - Tp^m- (33) 

The renormalization prescription for Ao specified by Eq. (10) is that YleAQ 1 {E) must vanish 
at M\ +2 — K 2 x /(2Mi2). This requires the inverse bare coupling constant to be 

= L (M 1+2 ) + ^k x . (34) 

Ao 27T 

Using this expression for 1/Ao and the expression for Lq(E) in Eq. (A13), the amplitude in 
Eq. (33) reduces to 



13 



+ 





+ 



FIG. 5: The connected Green's function iAo(E) for D\D 2 — > D\D 2 at th order in g can be 
obtained by summing a geometric series of one-loop diagrams. 

If we were to use the complex mass scheme, the amplitude Lq(E) would be given by 
Eq. (31) with M\ +2 replaced by M\ +2 — iT 2 /2. If we again demand that Re^Q 1 ^) must 
vanish at M\ +2 — k x /(2Mi 2 ), the renormalized expression for the amplitude would be 



ME) 



2n/M- 
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7(E) - Re7x' 
where j(E) is a complex energy variable defined by 



j(E) = V-2M 12 (£ - M 1+2 + iY 2 /2) 
and 7x is its value at E = M\ +2 — k x /(2Mx2): 



k\ - iM 12 T 2 . 



(36) 



(37) 



(38) 



The variable j(E) vanishes at the complex threshold E = Mi +2 — iT/2 and it has a negative 
imaginary part for all real values of E. 



B. Propagator for X 

If the local composite operator \qD\D 2 (x) is used as an interpolating field for X, the 
propagator for X is given in Eq. (24). The diagrams for the propagator of X at th order in 
g are shown in Fig. 6. These diagrams form a geometric series whose sum in the rest frame 
P = is 

i\ 2 L (E) 



A x (E,0) 



(39) 



1-X L (EY 

Using the expression for Aq(E) in Eq. (33), we can express the propagator for X as 

A x (E,0) = -t[X + A (E)}. (40) 

To obtain a renormalized propagator for X that does not depend on the ultraviolet cutoff, 
an additive renormalization is necessary. After adding the constant i\q, we obtain the 
renormalized propagator —iAo(E). 

The term — iAo(E) in the propagator for X in Eq. (40) has a pole at the complex energy 
Epoie given in Eq. (11). Using the renormalized expression for Aq(E) in Eq. (35), we find 
that the binding energy and width of X at th order in g are 

E x = 4/(2M 12 ), 
T x = 0. 



(41a) 
(41b) 
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FIG. 6: Feynman diagrams for the X propagator at th order in g. The interpolating field for the 
X is A £>i£> 2 - 



The wavefunction normalization factor for X is 

2tt 



M? 2 



kx. (42) 



If we were to use the complex mass scheme, the amplitude Aq(E) would be given by 
Eq. (36). This amplitude has a pole at a complex energy E po \ e that determines the binding 
energy and width of X. The complex parameter -E po i e satisfies 

7 (£ polc )-Re 7x = 0, (43) 

where j(E) and 7x are defined in Eqs. (37) and (38). Comparing with Eq. (11), we find 
that at LO in the complex mass scheme, the binding energy and width of X are 

E x = (Re 7 x)7(2M 12 ), (44a) 
T x = T 2 . (44b) 

If we treat T 2 as being of order g 2 , the binding energy of X in Eq. (44a) differs from 
k 2 x /(2M\2) only at order g 4 . The wavefunction normalization constant for X at LO in the 
complex mass scheme is 

2-7T 

Zx = ME Reix - (45) 



C. Short-distance Production of X 

The operator product expansion of the T-matrix element for a short-distance production 
process A — > B + X is given in Eq. (28a). The vacuum-to-X matrix element can be 
calculated via the LSZ formalism using X DiD 2 (x) as an interpolating field for X [24]: 

(X|A OJ D^(O)|0) = 4 /2 . (46) 

At th order in g, the normalization constant Zx is given in Eq. (42). At LO in the complex 
mass scheme, Zx is given in Eq. (45). The T-matrix element in Eq. (28a) can be expressed 
as the product of a short- distance factor and a long-distance factor: 

T[A -> B + X] = y / 2M x C B / 2 Z][ 2 . (47) 

The rate for producing X is obtained by squaring the T-matrix element in Eq. (47) and 
integrating over the appropriate phase space. If A consists of a single particle, its decay rate 
into B + X can be expressed in the factored form 

T[A^B + X] = T B A M 12 \Z X \. (48) 
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We have followed Ref. [24] in choosing the long-distance factor in Eq. (48) to be M\ 2 \Zx\, 
which is dimensionless. The short-distance factor is 

where the 4- vector P M satisfies P 2 = M?, ,. 



D. Short-distance Decay of X 

The operator product expansion of the T-matrix element for a short-distance decay pro- 
cess X — > C is given in Eq. (29). The X-to-vacuum matrix element can be calculated via 
the LSZ formalism using \ D\d\(x) as an interpolating field for X [24] : 3 

(®\\ D 1 D 2 (0)\X) = Z x /2 . (50) 

At th order in g, the normalization constant Zx is given in Eq. (42). At LO in the complex 
mass scheme, Zx is given by Eq. (45). The T-matrix element in Eq. (29) can be expressed 
as the product of a short- distance factor and a long-distance factor: 

T[X -> C] = yf2Mx~ Cg Z][ 2 . (51) 

The decay rate of X into the particles represented by C is obtained by squaring the 
T-matrix element in Eq. (51) and integrating over the phase space of those particles. It can 
be expressed in the factored form 

T[X -> C] = T c M 12 \Z X \. (52) 

We have followed Ref. [24] in choosing the long-distance factor to be Mi 2 \Zx\, which is 
dimensionless. The short-distance factor T c is 



T c = _L f TT ^Pi \ C g\ 2 (2n) A 5^(P - 5». (53) 



where the 4- vector satisfies P 2 = Mf +2 . 



E. Line shape of X in a Short-distance Decay Channel 

If X can be produced via the short-distance process A — > 5 + X and if it can decay via 
the short-distance process X — > C, then there can be resonant enhancement of the process 
A ^ B + C, where C represents the same particles but with a variable invariant mass 
M instead of Mx- The T-matrix element for this process can be expressed as the double 



3 Our notation might suggest that the matrix element in Eq. (50) is the complex conjugate of the matrix 
element in Eq. (46). However \X) in Eq. (50) is an in state, while (X\ in Eq. (46) is the hermitian 
conjugate of an out state. These two states are related by the S-matrix: \X, out) = S\X, in). 
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operator product expansion in Eq. (30). The Fourier transform of the matrix element is just 
the X propagator evaluated at (E,P) = (M, 0). The T-matrix element therefore reduces 
to 

n^B + C\=C? + mJ&gfa. (54) 

The T-matrix element can be expressed in a form in which the short-distance effects and 
long-distance effects are separated: 



T[A -> B + C] = -iC*> 12 Cg [Ao(M) - {2ir/M 12 )c^\ , (55) 

where Aq(E) is given in Eq. (35) and complex constant with dimensions of length 

that is completely determined by short-distance factors: 

b,c _ - M 12 {C^ C - ^ C^' 12 C^) 

Ca ~ ~ l • (56) 

The nonresonant term c^ ,C in Eq. (55) is a combination of short-distance factors, so the 

B C 

natural scale for c A ' is 1/m. The condition Kx <C m and the condition on M in Eq. (8) 
imply that the nonresonant term in Eq. (55) is small compared to the resonant term in the 
threhsold region secified by Eq. (8). We therefore set c^ ,c = 0. 

The invariant mass distribution of the particles in C is obtained by squaring the T-matrix 
element and integrating over the momenta of all the particles in the final state. If A consists 
of a single heavy particle, the invariant-mass distribution of the particles in C has the form 

1L[A^B + C\ = (r*r c ) H |^(M)| 2 . (57) 

The short- distance factors and r c are the same as in Eqs. (48) and (52). At th order 
in g, Ao(E) is given in Eq. (35). At LO in the complex mass scheme, Aq(E) is given in 
Eq. (36). Note that if the invariant mass distribution in Eq. (57) is divided by the product 
of the decay rates T[A — > B + X] and T[X — > C] in Eqs. (48) and (52), the short- distance 
factors cancel. This combination of observables has only a long-distance factor given by 
|A(M)| 2 /(27t|Z x | 2 ). 

The line shape of X as a function of the invariant mass E of the particles in C is 
given by the factor |^4 (i?)| 2 in Eq. (57). In Fig. 7, we compare the line shapes at th 
order in g and at LO in the complex mass scheme. At th order in g, the line shape 
diverges at Mi +2 — k, x /(2M 12 ). In the complex mass scheme, the line shape is a resonance 
whose maximum is near Mi +2 — k x /(2Mi 2 ) and whose width is determined by T 2 . If T 2 
k 2 x /(2M 12 ), the line shape for \E - M 1+2 + k 2 x /(2M 12 )\ < k x /(2M 12 ) is a nonrelativistic 
Breit-Wigner resonance whose full width at half- maximum is T 2 : 

\A ( E) I 2 ~ ** 2 "x/ M t2 f5g x 

lME)] ~[E-M 1+2 + «y(2M 12 )P + rl/4- (58) 

However, a Breit-Wigner resonance has tails that fall off like E~ 2 and it is integrable. In 
contrast, the line shape |^4o(-^)| 2 with Aq{E) given by Eq. (36) has tails that fall off only 
like E~ x and it is therefore not integrable. Thus the area under the resonance is not as 
well-defined as for a Breit-Wigner resonance. 



B,C 
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FIG. 7: Line shapes l^o^)! 2 of the X resonance in a short-distance decay channel at leading order 
in g. The dotted line is at order g° in ordinary perturbation theory (using Aq(E) in Eq. (35)). The 
dashed lines are at LO in the complex mass scheme (using Aq(E) in Eq. (36)). The line shapes 
are shown for k x = 27.8 MeV and for both T 2 = 0.07 MeV (upper solid line) and T 2 = 0.28 MeV 
(lower solid line). The units on the vertical axis are 10 -5 MeV -4 . 



IV. SECOND ORDER IN g 

In this section, we calculate the results in Section III to next-to-leading order in g. Since 
the D 2 — Di<p interaction allows transitions to D\D\§ states, we also calculate the rates for 
processes whose final state includes D\D\§. 



A. Amplitude for D 1 D 2 — > D 1 D 2 

One of the basic building blocks for the amplitudes for processes near the D 1 D 2 threshold 
is the amplitude for the propagation of E)\D 2 between contact interactions. The term in 
this amplitude of order g°, iLo(E), is given in Eq. (31). The term of order g 2 , iL 2 (E), is the 
sum of the three Feynman diagrams in Fig. 8. These diagrams are calculated in Section A 4 
of the Appendix. The final expression for L 2 (E) is the sum of the amplitudes L 2a (E) and 
L 2b (E) given by Eqs. (A21) and (A29): 

Ml 9 

L 2 (E) = L 2 (2M 1 + m) + _^F(k(£)), (59) 
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(a) 



(b) 



(c) 



FIG. 8: Diagrams of order g 2 for the amplitude for the propagation of D\D 2 between contact 
interactions. The open dots indicate that the vertex factors — zAq are omitted. 



where k{E) is the energy variable denned in Eq. (32), the function F(k) is 



F(k) 



g 2 M 12 m 1 



IT' 



1 dz 



mi 



tHz) 



M 



12 



l-t(z) 

1/2 



1/2 



In 



zk 2 + (1 - z)(m 1 /M 12 )(K 2 - k\) 



ZK7 



2 In 



2k 



2/?i k 



In 



K — «i 



and Ki is the value of k(E) at the D\Di(j) threshold 2M\ + m: 

k x = [2M 12 (M 2 _! -m)] 1 / 2 . 



(60) 



(61) 



The function t(z) in the integrand in Eq. (60) is a rational function of z that increases from 
to 1 as 2 increases from to 1: 

z 2 

= (2- z) 2 -4(m 1 /m) 2 (l- z) 2 ' ^ 

The term L 2 (2Mi +m) in Eq. (59) is real valued and includes logarithmic ultraviolet diver- 
gences. The function F(k) vanishes at the D\D\§ threshold k = K\ by definition. The first 
few terms in the Laurent expansion of F(k) around the E\D 2 threshold at k = can be 
obtained from Eqs. (A24) and (A31). The leading term is 



F(k) 



.g 2 M 2 2 {m 1 /M 12 fl 2 K 1 



2ttk 



(63) 



At th order in g, the ultraviolet divergence in Lq(E) was eliminated by summing a 
geometric series in Lq and then renormalizing the coupling constant Ao- The ultraviolet 
divergence in L 2 (E) can be eliminated in a similar way. The geometric series of diagrams 
that must be summed are those that can be obtained by replacing each one-loop subdiagram 
in Fig. 5 by the sum of the one-loop subdiagram and the three diagrams for L 2 (E) in Fig. 8. 
The sum of that geometric series gives an amplitude iA 2 (E) that can be obtained by making 
the substitution L — ► L + L 2 in Eq. (33): 



iA 2 {E) 



(64) 



1/X -[L (E) + L 2 (E)} 

Our renormalization prescription for Ao in Eq. (10) requires Re A 2 l (E) to vanish at M\ +2 — 
k|"/(2M 12 ). Using this prescription, the amplitude in Eq. (64) can be expressed as 



A 2 {E) 



2n/M 



12 



k(E) + F(k(E)) - Re [k x + F(k x - ie)\ 



(65) 
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FIG. 9: Line shapes |„4(i?)| 2 of the X resonance in a short-distance decay channel. The dotted lines 
are at order g 2 in ordinary perturbation theory (using A2(E) in Eq. (65)). The dashed lines are at 
LO in the complex mass scheme (using Aq(E) in Eq. (36)). The line shapes are shown for k,x = 27.8 
MeV and for both r 2 = 0.07 MeV (the two narrower resonances that are almost indistinguishable) 
and for r 2 = 0.28 MeV (the two wider resonances). The line shapes are normalized so their 
maximum values are 1. The inset shows their behavior near the D1D2 threshold. 



The quantity (^(-E")] 2 with A2(E) given in Eq. (65) is the line shape of the X resonance 
in a short-distance decay channel to order g 2 in ordinary perturbation theory. In Fig. 9, we 
compare this line shape with that at LO in the complex mass scheme, which is |^4o(-^)| 2 
with Ao(E) given by Eq. (36). They are shown for two values of the coupling constant 
g that correspond to the width of D 2 being T 2 = 0.07 MeV and 0.28 MeV. The two line 
shapes have the same qualitative behavior below the resonance and near the peak of the 
resonance. They have qualitatively different behavior near the DiD 2 threshold E = M 1+2 . 
Unlike |^4 (-^)| 2 ? the function |^4 2 (-E)| 2 vanishes at the DiD 2 threshold and has a sharp peak 
just above the threshold. The inset in Fig. 9 shows the behavior near the D 1 D 2 threshold, 
which is clearly unphysical. 

The reason l^^)! 2 vanishes at the DiD 2 threshold can be seen from the Laurent ex- 
pansion of F(k) around k = 0. The leading term, which is given in Eq. (63), diverges like 
1/k, as E — > Mi + 2- No matter how small the coupling constant g, the term proportional to 
g 2 1 'k, cannot be treated as a perturbation in the region near k = 0. The resolution of the 
problem can be found by noting that the diverging term in Eq. (63) can be expressed as 

— -<-^, (66) 

where T 2 is the width of the D 2 from its decay into D\4>, which is given in Eq. (13). Thus 
the divergence at k = implies that near the DiD 2 threshold, the imaginary part of the 
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.D 2 self-energy must be resummed to all orders. That resummation shifts the branch point 
in the leading order amplitude from the real threshold E = Mi +2 to the complex threshold 
E = M\ + 2 — 1^2/2. Making this change and then imposing our renormalization scheme, the 
amplitude in Eq. (65) becomes 

2vr/M 12 

2( } ~ 7 (E) + F(k(E)) + iM 12 T2/(2k) - Re[ 7 x + F(k x - ie)\ ' (b ° 

where j(E) and 7x are defined in Eqs. (37) and (38). This amplitude is a smooth function 
of E at k — 0. 

A convenient way to implement the resummation in Eq. (67) that can be extended staight- 
forwardly to higher orders in g is to use the complex mass scheme. This automatically gives 
an amplitude A.i{E) that has smooth behavior at the D1D2 threshold and is correct through 
relative order g 2 . In the complex mass scheme, the contribution to L 2 (-E) from the diagrams 
in Figs. 8(a) and 8(b) are obtained by replacing the energy variable k(E) by 7(E) defined 
in Eq. (37) and by replacing K\ by 7x, which is the value of j(E) at the D\D\§ threshold 
2Mi + m: 

71 = [2M 12 (M 2 _! -m- ^r 2 /2)] 1/2 . (68) 

In the complex mass scheme, there are additional contributions to the function L^E) coming 
from the one-loop diagram in Fig. 8(c) with the counterterm vertex in Eq. (21b). The 
additional contributions are given in Eq. (A33). Adding those terms is equivalent to making 
the substitution 

F(k) — F(k) + + (Z 2 ^ - 1) (^A,(M 1+2 ) + . (69) 

Using the expressions to order g 2 for T 2 in Eq. (13) and Z^ 1 in Eq. (20), we can see that 
the additional terms in Eq. (69) cancel the k~ x and k terms coming from the Laurent 
expansions of L 2 (,(£') in Eq. (A31). Our renormalization prescription for Ao in Eq. (10) 
requires ReA 2 ~ 1 (E) to have a zero at M i+2 — k 2 x /(2M 12 ). Using this prescription, the 
amplitude can be expressed as 

A ( F ) = ^l M ^ 2 ( 7(]) 

21 1 1 (E) + F cm ME))-^llx + F cms (lx)Y 1 1 

where j(E) and ^ x are given by Eqs. (37) and (38) and the function F cms (7) is 



F h) = g2Ml2mi 

- 1 cms V 7 } 2 



7T^ 



1 dz ( t\z) \ 1/2 / (l-^)(m 1 /M 12 )( 7l 2 -7 2 )-^ 2 _ ^ 
~ z 2 \l-t(z)J \ (l-z^mjMuh 2 



M l2 J V 7i 2 7i7 7i-7 



(71) 



The function F cms ('j) can be obtained from F(k) in Eq. (60) by subtracting several terms 
in the Laurent expansion in k and then replacing k(E) and K\ by ^y(E) and 71. The terms 
that are subtracted are the imaginary term of order /c -1 , the real terms of order k°, and 
the imaginary term of order k in the second of the two terms on the right side of Eq. (60). 
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FIG. 10: Line shapes |„4(i?)| 2 of the X resonance in a short-distance decay channel using the 
complex mass scheme. The dashed lines are at LO in the complex mass scheme (using Aq(E) in 
Eq. (36)). The solid lines are at NLO in the complex mass scheme (using A 2 (E) in Eq. (70)). 
The line shapes are shown for kx = 27.8 MeV and for both T2 = 0.07 MeV (the two narrower 
resonances that are almost indistiguishable) and for T2 = 0.28 MeV (the two wider resonances). 
The line shapes are normalized so their maximum values are 1. 



The real terms of order n° are subtracted because they would cancel in the denominator in 
Eq. (70) anyway. The leading terms in the Laurent expansion of Fans (7) around the complex 
D\D 2 threshold at 7 = are 



- 1 cms V I ) 



71 



/ \ 1/2 / AyT \ 12 

m . m l m l \ 1 M 12 \ 7 

— arcsm h + 2 ' 



mi m \Mi2 ) \mi J 71 



(72) 



For real values of E, there is an ambiguity in the function F cms (^(E)) in Eq. (70) for 
E < 2Mi + m because of the branch point at 7 = 71 in Eq. (71). This ambiguity can be 
resolved by the substitution 7 — > 7 — it. 

In Fig. 10, we compare the line shapes of X in a short-distance decay channel at NLO in 
the complex mass scheme, which is |^42(-E)| 2 with A 2 (E) given by Eq. (70), and that at LO in 
the complex mass scheme, which is | -4_o C^") 1 2 with Aq{E) given by Eq. (36). They are shown 
for two values of the coupling constant g that correspond to the width of D 2 being Y 2 = 0.07 
MeV and 0.28 MeV. The curves are normalized so their maximum values are 1. The line 
shapes at LO and NLO have the same qualitative behavior. As g decreases, the quantitative 
differences between the line shapes at LO and NLO decrease. Thus the complex mass scheme 
seems to provide a good solution to the infrared problem near the DiD 2 threshold that arises 
in ordinary perturbation theory. We will therefore use the complex mass scheme in most of 
the calculations in the remainder of this Section. 
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B. Propagator for X 

If the local composite operator \ Q DiD 2 (x) is used as the interpolating field for X, the 
propagator for X is defined in Eq. (24). The expression for the propagator at th in g is 
given in Eq. (39). A renormalizable expression for the X propagator that is accurate to 
order g 2 can be obtained simply by making the substitution L — > L Q + L 2 in Eq. (39): 

a ( p n) i# [L (E) + L 2 (E)] 

Ax{E ^-l-X [L (E) + L 2 (E)Y (73) 

Comparing with the expression for the amplitude A 2 (E) in Eq. (64), we see that the prop- 
agator in Eq. (73) can be expressed as 

A x (E,0) = -t[X + A 2 (E)]. (74) 

The expression for A 2 (E) at NLO in the complex mass scheme is given in Eq. (70). This 
amplitude has a pole at a complex energy E po \ c that determines the binding energy and the 
width of X. The energy -E po i e is the solution to 

7 (£ polc ) + F cms ( 7 (£ polc )) - Re [ lx + F cms ( 7x )] = 0, (75) 

where 7 (-E l ) and 7 x are defined in Eqs. (37) and (38). The wavefunction normalization 
constant for X is 

7 = 2n Re 7 (E pole ) 
X Mf 2 l + F c U 7 (£ polc ))- 1 > 

If we treat T 2 as order g 2 , the solution to Eq. (75) for the complex variable E po i e through 
order g 2 is 

E pole ^M 1+2 -K 2 x /(2M 12 )-i[r 2 /2-K X lmF cms (K X )/M 12 }. (77) 

Using the expression for -F cms ( 7 ) in Eq. (71), we find that the binding energy and width of 
the X to order g 2 are 

E x » 4/(2M 12 ), (78a) 



2g 2 m 1 n x 
1 2 H 



7T 



z 2 \l-t(z) J \M- 



'12 



(78b) 



where the upper endpoint zq of the integral is 



_ ( mi /M 12 )(K 2 -K 2 x ) 
Z - 7~ /nr \(..2 ..1 \ , ..2 • l' y J 



( mi /M 12 )(K 2 - K 



X) "T" K X 



For k x = 27.8 MeV, the second term on the right side of Eq. (78b) is 0.011 T 2 . Thus the 
width of the constituent D 2 accounts for most of the width of X at order g 2 . 

We can obtain an analytic expression for the width T x if k x (m 1 /M 12 ) 1 / 2 /ti. In this 
limit, the upper limit z on the integral approaches 1 and the width reduces to 



2g 2 m l K X 
1 2 H 



7T 



m .mi l vti\ ' 
arcsm 

1 "h 2 



mi m V M- 



(80) 



As m/Mi varies from to oo, the factor in square brackets ranges from it/ 2 to \ — ^M\ +2 /M 2 . 
Note that the expression for Tx in Eq. (80) can be smaller than T 2 for some values of the 
parameters. 
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C. Short-distance Production of X 

The operator product expansion of the T-matrix element for the short-distance produc- 
tion process A — > B + X is given in Eq. (28a). At NLO in the complex mass scheme, the 
operator matrix element is given by Eq. (46), where Zx is now the wavefunction normaliza- 
tion constant for X in Eq. (76). The factored expression for the T-matrix element is given 
by Eq. (47). If A consists of a single particle, its decay rate into B + X can be expressed in 
the factored form in Eq. (48), where Zx is given in Eq. (76). 



D. Short-distance decay of X 

The operator product expansion for the T-matrix element for the short-distance decay 
process X — > C is given in Eq. (29). At NLO in the complex mass scheme, the operator 
matrix element is given in Eq. (50), where Zx is the wavefunction normalization constant 
for X in Eq. (76). The factored expression for the T-matrix element is given by Eq. (51). 
The decay rate of X into the particles represented by C can be expressed in the factored 
form in Eq. (52), where Z x is given in Eq. (76). 



E. Line shape of X in a Short-distance Decay Channel 

The T-matrix element for the short-distance process A — > B + C, where C is a set 
of particles with invariant mass M close to the D 1 D 2 threshold, can be expressed as the 
double operator product expansion in Eq. (30). The T-matrix element at order g° is given 
in Eq. (54). An expression for this T-matrix element that is accurate through order g 2 can 
be obtained by substituting L — > L + L 2 : 

T[A^B + C]= C B A C + C^Cg (3D 
L J A A 12 1- X [L (M) + L 2 (M)} 1 ; 

The expression for the T-matrix element in which short-distance and long-distance factors 
are separated is given by Eq. (55) with Aq(M) replaced by A 2 (M). The invariant mass 
distribution for the short-distance decay channel is 

<T :[A ^B + C] = (T B A T C ) ^ \A 2 (M)\ 2 , (82) 



where and r c are the short-distance factors defined in Eqs. (49) and (53). The expression 
for A 2 (M) at NLO in the complex mass scheme is given in Eq. (70). The line shapes at LO 
and NLO in the complex mass scheme are compared in Fig. 10. 

If the invariant mass distribution in Eq. (82) is divided by the product of the decay rates 
F[A — > B + X] and T[X — > C] in Eqs. (48) and (52), the short-distance factors cancel, 
leaving only the long-distance factors \A 2 (M)\ 2 /(2ir\Z x \ 2 ), where A 2 (E) and Z x are given 
in Eqs. (70) and (76). The corresponding quantity at LO in the complex mass scheme is 
|„4, (^0| 2 /(27r|Zx| 2 )) where Ao(E) and Zx are given in Eqs. (36) and (45). The dependence 
of this long-distance factor on the invariant mass M is illustrated in Fig. 10, where the 
curves are all normalized to maximum value 1. The peak value of \A(M)\ 2 /(2tt\Zx\ 2 ) is 
also completely determined by long distances. For k x = 27.8 MeV and T 2 = 0.07 MeV, the 
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FIG. 11: Diagrams at order g for the vacuum-to-Di-D 2 (/> matrix element of the operator Ao-D{£> 2 . 



peak value decreases from 130.1 MeV" 2 at LO to 119.2 MeV" 2 at NLO. For k x = 27.8 MeV 
and T 2 = 0.28 MeV, the peak value decreases from 8.34 MeV" 2 at LO to 6.58 MeV -2 at 
NLO. The decrease in the difference between the peak values at NLO and LO from about 
27% at T 2 = 0.28 MeV to about 9% at T 2 = 0.07 MeV supports the assumption that the 
perturbative expansion is well-behaved in the complex mass scheme. 



F. Short-distance production of D\D\4> 

The fundamental theory may have short- distance production processes A — > B + D±Di(f), 
where A and B both represent one or more particles whose momenta in the D\D\<\) rest 
frame are of order m or larger. The D\D\(\) invariant mass M is assumed to be close to 
M 1+2 , as specified by the condition in Eq. (8). The operator product expansion of the 
T- matrix element for this process is given in Eq. (28b). 

We first calculate the rate for A — > B + D\D\(f) to order g 2 using ordinary perturbation 
theory. At order g, the operator matrix element in Eq. (28b) is given by the series of 
Feynman diagrams in Fig. 11 summed over the two permutations of the external Di lines. 
The sum of the geometric series of diagrams is Aq(M) / (— Ao) multiplied by the first diagram 
in the series. According to the Optical Theorem, the square of the first diagram integrated 
over the 3-body phase space for DiD\(j) is given by the sum of the 3-body cuts through the 
Feynman diagrams for Li2(E) in Fig. 8. If M < M1+2, the sum of the 3-body cuts gives the 
imaginary part of L 2 (E). Thus the square of the T- matrix element in Eq. (28b) integrated 
over the 3-body phase space of DiDi<p is 

= |Cj' 12 | 2 |A(M)| 2 (-2)Imi 2 (M), (83) 

where pi, p 2 , and k are 4- momenta for Dx, Di, and <f), respectively, and P 2 = M 2 . If A 
consists of a single particle, the invariant mass distribution for A — > B + DiD\(f) at order g 2 
can be expressed as 

^-[A ^B + DiDi4>] = T B A ^L \ MM )\ 2 (-2) ImL 2 (M), (84) 

where is the short-distance factor defined in Eq. (49) and we have replaced a multiplica- 
tive factor of M by Mi +2 . 

If there is a short-distance process A — » B + DiDi(f>, then the process A —>■ B + D\D 2 
is also allowed. The invariant mass distribution for A — ► B + DiD 2 at order g° is given 
by the expression on the right side of Eq. (84) with Im L 2 (M) replaced by ImL (M). The 
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threshold for this process is M > M 1+2 . Since D 2 ultimately decays into Di<f>, a decay into 
D\D 2 can also be regarded as a contribution to the inclusive decay into D\D\(j). 

The result in Eq. (84) is an order g 2 correction to the inclusive invariant mass distribution 
for D 1 Di<j). It is nonzero for all M > 1M\ + to. For M > M 1+2 , there are additional 
corrections of order g 2 that can be obtained by replacing any of the one-loop subdiagrams in 
Fig. 11 by one of the two- loop diagrams for L 2 (E) in Fig. 8. Nonperturbative renormalization 
of the coupling constant A then requires that a geometric series of these corrections be 
summed to all orders. The net effect is to replace Ao(M) in Eq. (84) by A 2 (M). Thus the 
complete result for the inclusive invariant mass distribution for D\D\<$> through order g 2 is 

<W [a^b + DM] = T B A ^\A 2 {M)\ 2 (-2)Im[L (M) + L 2 (M)], (85) 



where the functions A 2 (M), L (M), and L 2 (M) are given in Eqs. (65), (31), and (59). Note 
that taking the imaginary parts of L (M) and L 2 (M) eliminates the additive ultraviolet- 
divergent terms. A more compact expression for the distribution in Eq. (85) is 

<K [A^B + D 1 D 1 0] = Y B A ^ lmA 2 (M). (86) 



dM A 7T 

The inclusive invariant mass distribution in Eq. (85) when calculated using ordinary per- 
turbation theory has unphysical behavior at the threshold M = M 1+2 . The factor |.4. 2 (M)| 2 
has a zero at M — M 1+2 and a sharp peak just above that threshold as illustrated in Fig. 9. 
The term — 2ImL 2 (M) also has unphysical behavior at this threshold. The imaginary part 
of L 2 (E) in the region 2M 1 + to < E < M i+2 is given by the sum of Eqs. (A22) and (A30). 
The resulting expression for — 2lmL 2 (E) is 



7T Z 



9 Zo dz ( t\z) \ 1/2 + ( m l \ 1/2 («x - /, ( E) Y 



z 2 \l-t(z)J \M 12 J 2K lf i(E) 
where the upper limit z of the integral is 



37) 



= { mi /M 12 ){K 2 -K{E) 2 ) 
(mi/M 12 )(«?-«(S) 2 )+«(S) 2 " 1 ' 

The function in Eq. (87) diverges like 1/k(E) as E — > M\ +2 . As a consequence, near the 
D\D 2 threshold, corrections of higher order in g 2 that come from the D 2 self energy are 
not suppressed. This is the same problem that prompted us to introduce the complex mass 
scheme in Section IV A. 

We next consider the inclusive rate for A — > B + D\D\(j> in the complex mass scheme. At 
LO in the complex mass scheme, the invariant mass distribution for DiDi<p can be obtained 
from Eq. (86) by replacing A 2 (M) by Ao(M) in Eq. (36). This can be written more explicitly 
as 

where 'y(E) and 7x are defined in Eqs. (37) and (38). This expression has a nonzero 
imaginary part at all energies E, even below the D^D^ theshold. For E < 2Mi + to, the 
imaginary part is small and it is cancelled by higher orders in the complex mass scheme. At 
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FIG. 12: Line shapes of the X resonance at NLO in the complex mass scheme. The solid lines 
are for a short-distance decay channel ([^(-E 1 )) 2 using A2{E) in Eq. (70)). The dash-dotted lines 
are for the D\D\^> channel (ImA2(E) using Ai{E) in Eq. (70)). The line shapes are shown for 
kx = 27.8 MeV and for both T2 = 0.07 MeV (narrower resonances) and T2 = 0.28 MeV (wider 
resonances). The line shapes are normalized so their maximum values are 1. 

NLO in the complex mass scheme, the invariant mass distribution for DiDi<fi is obtained 
from Eq. (86) by using the expression for A 2 (M) in Eq. (70). 

In Fig. 12, we compare the line shape of X in the D\D^<\) channel, which is proportional to 
lmA2{E), with the line shape of X in a short-distance decay channel, which is proportional 
to |„4 2 (-^)| 2 - They are shown for two values of the coupling constant g that correspond to 
the width of D 2 being T 2 = 0.07 MeV and 0.28 MeV. The curves are normalized so their 
maximum values are 1. The line shapes have the same qualitative behavior. The line shape 
for D\Di(f) is shifted towards higher energy by an amount that decreases as g decreases. 

The Belle Collaboration has observed an enhancement in the production of D D tt near 
threshold in the decay B — > K + D°D°-k° [32]. The enhancement peaks at the invariant 
mass 3875.2l};g MeV, where we have combined the errors in quadrature. If we use the recent 
precision measurement of the D° mass by the CLEO Collaboration [26], the peak is above 
the D*°D° threshold by 3.4^'g MeV. In contrast, the peak observed in the short-distance 
decay channel <//■?/> 7r + 7r~ is below the D*°D° threshold by 0.6 ± 0.6 MeV. The difference 
is 4.0^2:0 MeV, which differs from by two standard deviations. Part of the discrepancy 
may be due to a shift in the invariant mass distribution analogous to the one illustrated in 
Fig. 12. 
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G. Decay of X into D X D X 4> 



The width Fx of the X resonance was defined in Eq. (11) in terms of the energy E vo \ c of 
the pole in the X propagator: Y x = — 2Im_E po i c . At LO in the complex mass scheme, Y x 
is simply T 2 . At NLO in the complex mass scheme, Yx is obtained by solving Eq. (75) for 
.Epoie- If the width of the X resonance is sufficiently small, its width can also be interpreted 
as its decay rate. A direct calculation of that decay rate will give a result that is closely 
related to but not identical to Y x . We will denote the result of the direct calculation of the 
inclusive decay rate by Y[X — > DiDi4>}. 

We can use the Optical Theorem to calculate the inclusive decay rate for X — > DiDi4>: 

Y[X - D 1 D 1 c/ ) } = ^-ImT[I - X}. (90) 

The forward T-matrix element for X — > X can be defined using the LSZ formalism and it 
is given in Eq. (27). At NLO in the complex mass scheme, the normalization constant Z x 
and the propagator A x (E,0) are given in Eqs. (76) and (73). The resulting expression for 
the decay rate of X in Eq. (90) is 

Y[X -> D 1 D 1 <f>] = -21m r Zx ^ Mx ] (91) 

We have used Eq. (64) to express this in a form with a factor of A 2 1 (M X ), because A 2 (E) 
does not depend on the ultraviolet cutoff A. The decay rate in Eq. (91) depends on A through 
the bare coupling constant A and through the ultraviolet-divergent additive constants in 
Lq(M x ) and L2(M X ). From the renormalized expression for A 2 (E) in Eq. (64), we can see 
that the divergent terms in L (E) + L 2 (E) must be cancelled by 1/Ao as A — > oo. The 
leading divergence is a linear divergence in L (Mi +2 ) that is linear in A. Thus 1/Ao must 
include a cancelling linear divergence. This implies that X [L (E) + L 2 (E)] must approach 
1 as A — > oo. Upon taking the limit A — > oo, Eq. (91) reduces to 

Y[X -> D 1 D 1 (f)} = -2 Im (Z X A 2 \M X )) . (92) 

An explicit expression for the inclusive decay rate of X can be obtained from Eq. (92) 
by inserting the renormalized expression for A 2 (E) in Eq. (70): 

Y[X - D 1 D 1 c^] = ^ (-2) Im [Z x (>y(M x ) + F cms ( 7 (M x )) - Re[ 7 x + E cms (lx)})} , 

(93) 

where Z x , j(E), r y x , and -F cms (7) are given in Eqs. (76), (37), (38), and (71). If we treat Y 2 
as order g 2 , the decay rate in Eq. (93) reduces at order g 2 to the corresponding result for 
the width Y x in Eq. (78b). 

In Fig. 13, we compare the width Y x obtained by solving Eq. (75) and the decay rate 
Y[X — > DiDi4>] in Eq. (93) as functions of Y 2 . They agree for small Y 2 , because the difference 
between Y x and Y[X -> DiDi<f>] is of order g 4 . They differ by less than 1% if Y 2 < 0.16 
MeV. 
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FIG. 13: The width Tx of the X resonance (solid line) and the decay rate T[X — > DiDi4>] (dashed 
line) at NLO in the complex mass scheme as functions of F 2 . The vertical axis is in units of T 2 . 
The intercept on the vertical axis is approximately 1.011. 

V. SUMMARY 

We have examined the effects of a weakly-bound hadronic molecule on a nearby 3-body 
threshold, where the three bodies consist of a primary constituent of the molecule and the 
decay products of the other constituent. We studied these effects in the simplest possible 
model: the scalar meson model, with spin-0 constituents Di, D 2 , and <fi and momentum- 
independent interactions. The primary constitutents of the molecule X are D\D 2 and the 
3-body threshold is that for DiDi<p. The D\D 2 contact interaction with coupling constant 
Ao was treated nonperturbatively. The D 2 — Di<p interaction with coupling constant g was 
treated perturbatively. Several observables were calculated to all orders in A and to next- 
to-leading order in g. Both ultraviolet problems and infrared problems were encountered in 
these calculations. 

We found that at next-to- leading order in g, the ultraviolet divergences can be removed 
by a perturbative renormalization of the D 2 mass, a nonperturbative renormalization of Ao, 
and a nonperturbative renormalization of short-distance coefficients in the operator product 
expansion. The nonperturbative renormalizations required the summing of a geometric series 
of order-g 2 subdiagrams to all orders. The subdiagrams represent order-g 2 contributions 
to the amplitude for the propagation of DiD 2 between contact interactions. Thus the 
renormalized amplitudes at next-to-leading order include all terms of order g 2 together with 
subsets of terms from all higher orders in g. 

We found that the next-to-leading order amplitudes suffer from an infrared problem at 
the D\D 2 threshold. The problem is illustrated dramatically in Fig. 9, which shows the 
line shape |„4(-E)| 2 in a short-distance decay channel as a function of the energy E. The 
pathological behavior at E — M 1+2 arises because the D 2 — Di<p interaction shifts the D 1 D 2 
threshold from the real value M\ +2 to the complex value M 1+2 — iT 2 /2. Since T 2 is of order 
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g 2 , a strict expansion in powers of g leads to singularities at E — M i+2 . One possible solution 
to this problem is to sum the geometric series of self-energy corrections to the D 2 propagator 
to all orders, but this leads to very complicated loop integrals. A simpler solution is to use 
the complex mass scheme, in which the width T 2 is included in the Feynman rule for the 
D 2 propagator and its effects are systematically compensated at higher orders in g through 
counterterms. This partial resummation of the D 2 propagator corrections leads to smooth 
behavior of the line shape |^4(i?)| 2 at the D Y D 2 threshold as illustrated in Fig. 10. Having 
solved the infrared problems, we calculated several observables to NLO in the complex mass 
scheme. They include the line shapes of the X resonance in both a short- distance decay 
channel and in the D\Di(j) channel. 

The X(3872) is a weakly-bound hadronic molecule whose primary constitutents are a 
superposition of charm mesons: D*°D° + D°D*°. Its mass is near the 3-body thresholds for 
D°D°ir , D°D~ir + , and D + D°ir~. The methods we applied to the scalar meson model can 
be extended straightforwardly to this multichannel problem. A conventional perturbation 
expansion in the D* — Dir coupling constant will have infrared problems at the D*°D° 
threshold which are associated with the decay D*° — > D 7t°. These infrared problems can 
be avoided by using the complex mass scheme. The ultraviolet problems can be expected 
to be more severe in this system, because the D* — Dn interaction is proportional to the 
3-momentum of the pion. Thus renormalization may not be as simple as in the scalar meson 
model. If the renormalization problem can be solved, it will be straightforward to calculate 
observables to NLO in the complex mass scheme. One of the most interesting applications 
will be to predict the difference between the line shapes of the X(3872) resonance in a 
short- distance decay channel, such as J/ip 7r + 7r~, and in the D D°ir channel. 
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APPENDIX A: LOOP INTEGRALS 

In this appendix, we calculate the loop integrals that appear at order g° and at order g 2 
in the scalar meson model. 

1. Dimensional regularization 

Many of the loop integrals we need to evaluate are ultraviolet divergent. They can be 
regularized by using dimensional regularization of the integrals over the 3-momenta. We 
analytically continue the 3-dimensional integrals to 3 — 2e dimensions using the prescription 




(Al) 
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where \i is the renormalization scale and (z) a is the Pochhammer symbol: 



The function of e in the prefactor is designed to simplify the analytic expressions for dimen- 
sionally regularized integrals by cancelling the effects of the analytic continuation of angular 

— » 

integrals. Thus the integral of a scalar function of k is 

jf/(|fc|)=Bm ^/i 2e J™dkk 2 - 2 J(k). (A3) 

Loop integrals are evaluated by first inserting the appropriate nonrelativistic propagators: 

i 

Po - Mi -p 2 /(2Mi) +ie 

where i = 1, 2, and = m. Integrals over the loop energy are then evaluated using the 
residue theorem. After using the Feynman parameter trick to combine denominators, the 
dimensionally regularized integrals over the 3-momenta can be evaluated. The final steps are 
the evaluation of the Feynman parameter integrals and the analytic continuation to e = 0. 
It will be convenient to express the results in terms of the variables 

k = [-2M 12 (E - M 1+2 ) - is} 1 ' 2 , (A5a) 
«i = [2M 12 (M 2 _! -m)] 1/2 . (A5b) 



2. Self-energy of D 2 



The one-loop diagram for the self-energy E(po,p) of D 2 is shown in Fig. 2(a). The 
expression for the diagram is 

f d 4 q 

Z 2 a(Po,p) = ~ig 2 J Ai(g , q) A<f(Po ~ Qo, W~ q\)- ( a 6) 



The residue theorem can be used to evaluate the energy integral: 

1 



S 2 a(P0,P) 



9 Lpo -M 1 -m- 



q 2 /(2M 1 ) - (p - q) 2 /(2m) + is' 



(AT) 



If we implement dimensional regularization using the prescription in Eq. (Al), the analytic 
result is 



g 2 m 1 /j, 2e 
27rcos(7re) 



-2m! 



Po — Mi — m 



2(M l + m) 



is 



(A8) 



The self-energy in dimensional regularization is given by the analytic continuation of 
Eq. (A8) to e = 0: 



^2a{po,p) 



g 2 mi 
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—2m\ 
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(A9) 
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The momentum integral in Eq. (A7) has a linear ultraviolet divergence that is set to zero 
by dimensional regularization. The self-energy in a general regularization scheme is given in 
Eq. (16). The linear ultraviolet divergence is contained in the extra term S 2(l (Mi + to, 0), 
which is real valued. 



3. Propagation between contact interactions at leading order 

The amplitude for the propagation of DiD 2 between contact interactions is given at 
order g° by the one-loop diagram in Fig. 4. In the center-of-momentum frame where the 
4-momentum is (E,0), this amplitude is 

/j4 
^LA 2 (q ,q)A 1 (E-q ,q). (AlO) 

The residue theorem can be used to evaluate the energy integral: 

LoiE) = l(E-M 1+2 )- q y2M 12 + i e- (AU) 

If we implement dimensional regularization using the prescription in Eq. (Al), the analytic 
result is 

Lo(E) = ^r-. 1 ^ (A12) 
27rcos(e7r) 

where k is defined in Eq. (A5a). The amplitude in dimensional regularization is given by 
analytic continuation to e = 0: 

L (E) = ¥* K . (A13) 

The momentum integral in Eq. (All) has a linear ultraviolet divergence that is set to zero 
by dimensional regularization. The amplitude in a general regularization scheme is given in 
Eq. (31). The linear ultraviolet divergence is contained in the extra term L (M 1+2 ), which 
is real valued. 



4. Propagation between contact interactions at next-to-leading order 

The amplitude for the propagation of D 1 D 2 between contact interactions has contribution 
of order g 2 from the two- loop diagrams in Fig. 8. In the center-of-momentum frame where 
the 4-momentum is (E,0), the order-g 2 amplitude iL 2 (E) is a function of E only. We 
proceed to calculate the contributions from each of the three diagrams. 



a. Correction from 4> exchange 



The contribution to the amplitude iL 2 {E) from the two- loop diagram in Fig. 8(a) in 
which is exchanged between the D 1 and D 2 is 

iL 2a (E) = -g 2 J-0- 4 J ^A 2 (k ,k)A 1 (E-k ,k) 

xA 2 (? , q) ME - g , q) A^(k + q -E,\k + q\). (A14) 
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The integrals over the loop energies ko and go can be evaluated by closing the contours in 
the upper half-planes. After introducing Feynman parameters and evaluating the integrals 
over the loop momenta, the amplitude reduces to 

L2a{E) = — pr- 1 » 

x [ dx [ X dy[(l-x)(l-y)-(m l /m) 2 (l-x-y) 2 ]- 3/2+e 
Jo Jo 

x [(x + y)K 2 + (l-x- y)( mi /M 12 )(K 2 - k\)] ~ 2e , (A15) 

where k and K\ are defined in Eqs. (A5). One of the Feynman parameter integrals can be 
evaluated analytically by changing variables to z = x + y and w = (x — y)/(x + y). The 
integral over w gives a hypergeometric function: 

2tt e J 

x 2 F 1 (| - 6, i, |; t(z)) [zn 2 + (1 - z)(m 1 /M 12 )(« 2 - k 2 )} ~ 2e , (A16) 

where t(z) is the rational function of z given in Eq. (62). The pole at e = in Eq. (A16) is 
a logarithmic ultraviolet divergence. 

The amplitude L 2a (E) in dimensional regularization is defined by the Laurent expansion 
of the expression in Eq. (A16) to order e°. The expansion in powers of e is facilitated by 
using a transformation formula to replace the hypergeometric function in Eq. (A16) by one 
of the form 2 .Fi(l, e, §! t)'- 

2^(1 - 6, i, |; f) = 2^(1, e, |; f)(l - f)" 1/2+e . (A17) 

The expansion of 2 -^i(l, e, |; to order e is simple: 

2^(1, e, §; t) = 1 + let 2 F 1 (1, 1, §;*) -h (A18) 

Expanding the amplitude L 2a (-E') in Eq. (A16) to th order in e, it reduces to 

L 2a £ = - 12 -arcsin^ --4 + 7 + 21n2 
27r d \ 2 m \ e 



+ mi Z" 1 cte ( t 3 (z) 



m J z 2 \ 1 — 

-2 In 

We have used the integral 

*dz ( t 3 (z) ^ 1/2 



t(z) 3 
^ 2 + (l-^)(m 1 /M 12 )(/€ 2 -/€ 2 ) 



/i 2 



(A19) 



o z 2 \1 — t(z) J 2m\ m 



m . mi . . 

arcsm — . (A20) 
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This can be evaluated by changing variables to w = z/[z + 2r(l — z)\. Most of the terms in 
Eq. (A19) are independent of E. The expression can be simplified by expressing it in the 
form 

T ,^ T x g 2 M 2 m x f 1 dz ( t 3 (z) N 1/2 

L 2a (£) = L 2a (2M 1 +m) + IJ ' ' 



7T 3 7 z 2 Vl ~t(z) 

yAn zn 2 + {I - z){ mi /M 12 ){K 2 - k\) ^ 
zk\ 

The amplitude L 2a {2M\ + m) at the D\D\<f) threshold is real- valued and includes the log- 
arithmic divergence. The amplitude L 2a (E) has an imaginary part for E > 2M X + m. For 
energies in the range 2M\ + m < E < M 1+2 , the imaginary part is 

g 2 M^m 1 r ZQ dz_ f t\z) ^ 1/2 
1^~J *\T=tCz) 



lmL 2a (E) = -^f^ 5(t^) , (A22) 



where the upper limit z of the integral is 

20 " (m 1 /M 12 )( K 2 -^ )+/t 2- ^) 
The first few terms in the Laurent expansion of L 2a (E) in powers of k are 

1/2 



L 2a (£)^L 2a (2M 1+ m) + ^i 2 ™ 1 



2tt 3 



t(*) y ' zm : 



12 



m .mi . / Mi2 \ k 
—t7r — arcsm h 2m 



mi m \ m i J K \ 



(A24) 



b. Correction from Di self- energy 

The contribution to iL 2 (E) from the two-loop diagram in Fig. 8(b) with a self-energy 
correction to the D 2 line is 

iL 2h {E) = -i J -0^ A 2 (k , kf A 1 (E - ho, k) £ 2a (A; , k), (A25) 

where E 2a (/co,/c) is the contribution to the D 2 self-energy from the one-loop diagram in 
Fig. 2(a). The analytic result for £ 2 a(&o> k) in dimensional regularization is given in Eq. (A8). 
Since £ 2a (/c , k) has a branch cut for ko in the lower half-plane, it is convenient to close the 
k contour in the upper half-plane. After introducing Feynman parameters, the integral over 

— * 

k can be evaluated. If we implement dimensional regularization using the prescription in 
Eq. (Al), the result is 

g'M^m, f mi \ 1/2 - e (l-2e) 2 (l) 2e 4e 



L2b(E) ~ 16tt3 \M 11 ) (i) 2 cos 2 (7re)e /i 



x 



fdxx(l - x)-^ 2+€ [xk 2 + (1 - x){M u /M 12 ){k 2 - K\)Y 2t , (A26) 
Jo 
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where k and K\ are defined in Eqs. (A5) and Mn is given by 

M 1 (M 1 + m) 



Mn = 



2Mi +m 



(A27) 



The pole at e = in Eq. (A26) is a logarithmic ultraviolet divergence. 

The amplitude L 2b (E) in dimensional regularization is defined by the Laurent expansion 
of the expression in Eq. (A26) to order e°. Expanding to th order in e, the expression 
reduces to 



L 2b (E) = 



g M( 2 mi ( nil 



47T 3 



M 



11 



1/2 



I-4 + 41n2-ln^- 



21n^ 



+J / 1 fe x(l-x)-^ 1 . 1 - 2 + ' 1 -" M "^" K2 - K " 
Jo K 



(A28) 



The integral over x can be evaluated analytically. We can simplify the expression by replacing 
M11/M12 by 1, because the difference is suppressed by a factor of k\/M 2 . Most of the terms 
in Eq. (A28) are independent of E. The expression can be simplified by expressing it in the 
form 



L 2b {E) = L 2h (2M 1 + m)- 



g Mf 2 mi f mi 



2tt 3 



M 



12 



!/ 2 / / \2 

2 In - 1 + ( K ~ Kl > \ n K + Kl 



2ki 



2kki 



K — Ki 



(A29) 



The amplitude L 2 t,(2Mi + m) at the DiDi<f) threshold is real-valued and includes the log- 
arithmic divergence. The amplitude L 2 b(E) has an imaginary part for E > 2Mi + m. For 
energies in the range 2Mi + m < E < Mi +2 , the imaginary part of L 2 b(E) is 



lmL 2b (E) 



g 2 M\ 2 mi ( mi \ 1/2 (k x - k) 2 



4vr 2 



M- 



12 



KKi 



(A30) 



Note that this diverges at the E X D 2 threshold k = 0. The first few terms in the Laurent 
expansion of L 2a (E) in powers of k are 



L 2b (E) w L 2b {2Mi + m) 



fM 12 mi / m 



2tt 3 



M 



12 



1/2 



1 - 2 In 2 + iir 



2k 



1 + 



K 

2~iTi 



(A31) 



c. Correction from counterterms 



In a general regularization scheme, the one-loop self-energy subdiagram S2 a (^o,^) i n 
Eq. (A25) has a linear ultraviolet divergence, which is included in the term Yt 2a (Mi+m, 0) in 
Eq. (16). The resulting divergence in L 2b (E) is cancelled by the one-loop diagram with a D 2 
mass counterterm in Fig. 8(c). The counterterm vertex is —i5M, with 5M = H 2a (M 1 +m, 0). 

In the complex mass scheme, the D 2 propagator counterterm has the more complicated 
form given in Eq. (23). With this counterterm, the contribution to L 2 (E) of order g 2 from 
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the one-loop diagram in Fig. 8(c) is 

iL 2c (E) = -i J ■0^A 1 (E-k o ,k)A 2 (k o ,k) 2 ((5M + tT 2 /2) 

+ {Z 2 l - l)[k - M 2 - A; 2 /(2M 2 )]). (A32) 

In the term with the factor Z^ 1 — 1, the integral is proportional to Lq(E). This amplitude in 
a general regularization scheme is given in Eq. (31). In the term with the factor 5M + iT 2 /2, 
the integral is convergent. The complete result for the diagram is 

L 2c {E) = (5M + *r 2 /2)|| + (Z, 1 - 1) (L (M 1+2 ) + . (A33) 
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